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Proposed  Objectives  of  the  Research 

The  objective  of  the  proposed  research  was  to  expand  ^  new 
semiempirical  method,  SAMI,  that  explicitly  includes  treatment  of  ^-o^dab. 
The  elements  for  which  parameters  were  to  have  been  deveioped  are  sulfur 
(S),  silicon  (Si),  phosphorous  (P),  arsenic  (As),  and  copper  (Lu). 

Attained  Objectives  of  the  Research 

Elemental  parameter  sets  have  been  completely 
validated  for  sulfui  (S),  silicon  (Si),  and  phosphorous  f 

been  completed  on  SAMI  parameters  for  iron  (Fe),  although  final  results  are 
not  available.  Dr.  George  Famini,  the  technical  officer  for  the  project 
“Rested  that  we  suspend  work  on  arsenic  (As)  and  copper  (Cu)  m  favor  of 
iron,  a  request  with  which  we  complied. 

Statement  of  Work  Completed 

1.  M-olecular  Basis  S.et  for  Rarameterization  (MBSP)  updated  and 
expanded  for  silicon. 

2.  MBSP  updated  and  expanded  for  sulfur. 

i.  MBSP  updated  and  expanded  for  phosphorous. 

4.  MBSP  created  for  iron. 

5.  MBSP  created  for  copper. 

6.  SAMI  parameters  completed  for  silicon. 

7.  SAMI  parameters  completed  for  sulfur. 

8.  SAMI  parameters  completed  for  phosphorous. 

9  Recoding  of  AMP  AC  and  the  parameterization  program  to  allow  full 
and  correct  communication  and  interchange  of  data  between  the  two 

10.  Exploration  of  mulitiplicity  and  electronic  state  problems  surrounding 
the  description  of  the  iron  systems  in  theMBSP. 

11.  Derivation  of  a  stable  and  predictable  Configuration  Interaction  (C  ) 

process  for  use  in  the  parameterization.  i 

12.  Preliminary  grid  searches  completed  for  one-center /one-electron 

13  Preparation  of  data  and  presentation  of  results  at  the  American 

Chemical  Society  National  Meeting  in  San  Diego. 

Results  and  Discussion 

As  noted  in  Appendkes  1,  2,  ^d  3.  Highly 
parameters  for  silicon,  sulfur  and  phosphorous  have  been  comp  ®  , 

results  show  a  substantial  Improvement  over  preceding  ‘°"f. 

semiempirical  methods,  especially  for  hypervalent  systems  d-orbila  s 

included  with  SAMI  have  made  an  important  difference  m  this  context.  This 


advance  is  a  substantial  addition  to  the  capabilities  of  computational 
chemists,  allowing  them  to  study  many  industrially  imjportant  processes  that 
have  previously  been  treatable  only  by  high  level  ab  initio  approaches. 

Appendix  4  lists  preliminary  results  for  iron  within  the  SAMI 
protocol.  As  is  readily  evident,  these  results  are  not  up  to  the  par  expected  for 
methods  of  this  complexity  and  is  certainly  not  in  keeping  with  historical 
performance.  As  noted,  however,  these  results  are  preliminary.  They  are 
presented  in  this  context  to  illustrate  that  the  method  is  a  viable  approach  to  the 
problem  of  modeling  transition  metals  with  semiempirical  methodology.  It 
should  be  emphasized  at  the  outset  of  this  discussion  that  the  work  on  iron 
and  other  transition  metals  is  a  completely  new  area  using  Dewar-style 
semiempirical  theory.  Previous  experience  has  largely  been  with  closed-shell 
systems  (singlet  spin  state)  and  with  much  more  "regular"  chemistry.  This  is 
certainly  not  the  case  with  transition  metals.  Semiempirical  methods  found 
their  first  real  application  in  the  context  of  organic  molecules  and  reactions 
and  have  only  recently  moved  to  heavier  main  group  elements. 

The  results  presented  here  for  iron  are  for  a  very  simplistic  closed- 
shell,  restricted  Hartree-Fock  treatment  of  the  wavefunction.  We  began  with 
this  approach  (which  we  certainly  do  not  intend  to  finish  the  project  using)  as 
a  way  of  exploring  the  problems  that  will  be  encountered  in  the 
parameterization  of  elements  with  high  spin  states  requiring  extensive  use  of 
Configuration  Interaction  (Cl).  In  this  light,  our  efforts  have  been  very 
instructive  and  valuable.  A  number  of  obstacles  have  been  identified  and 
most  have  already  been  overcome.  Some  of  these  are  listed  below.  An 
asterisk  (*)  indicates  that  we  have  solved  the  problem  and  are  moving 
forward. 


•  re-coding  of  AMPAC  to  handle  very  large  open-shell 
multiplicities  such  as  are  routinely  encountered  with 
transition  metals  * 

•  development  of  new  Cl  algorithms  and  approaches  * 

•  errors  in  the  structure  of  the  parameterization  program  that 
led  to  incorrect  passing  of  data  between  it  and  AMPAC  when 
Cl  was  invoked  * 

•  difficulty  extrapolating  initial  values  for  the  parameters 
from  which  to  begin  searches  * 

•  energy  separation  of  the  parameters  describing  the  one- 
electron/one-center  energies 

•  ordering  or  the  atomic  orbital  levels 

•  partially  automated  approach  to  Cl  for  easier  calculation  and 
determination  of  the  minimum  level  of  Cl  required  for  a 
particular  system 


♦  lack  of  high  quality  experimental  data  on  spin  multiplicity  * 

The  implementation  of  "automatic  Cl"  will  be  an  important  consideration  in 
future  work  with  transition  metal  systems.  Calculations  on  organic  systems 
require  Cl  to  be  used  only  in  very  special  circumstances,  whereas  it  mU  nee 
lole  applied  almost  routinely  with  transition  metals.  We  intend  to  add  a  set 
of  program  code  that  will  predict  which  MOs  occupied  by  which  electrons 
3  to  be  mixed  and  automatically  provide  that  pattern  “amg  varioi^s 
MIPAC  options  (the  OPEN  keyword).  (The  user  will  obviously  be  able  to 
override  these  default  values  if  the  need  arises.)  At 

the  "natural"  first  guess  from  AMPAC  to  commence  the  SCF  procetoe  for 
automatically  assigning  the  Cl  mixing  pattern.  In  this  abroach,  toe  RH 
wavefunction  from  which  the  microstates  are  “"f 
mixing  is  one  where  orbital  occupancy  and  degeneracy  P® 
dSd  cTen-shell  systems  are  trl^ated  in  this  context  using  the  hah-electron 
(more  apt^named  the  "partial  electron")  approximation  pioneered  by  Dewar 

There  S' a  problem  with  this  from  the  parameterization  standpoint  in 
that  the  natural  pattern  first  generated  is  almost  exclusively  dependerit  on  the 
Ucc  and  Ddd  parameters.  This  means  that  the  "correct  answer  in  terms 

of  orbhal  degene^rLy  and  ordering  must  be  an  input  to  the  parameterization 
run  rather  than  a  result.  We  plan  to  overcome  this  limitation  by  applying 
SLr  S  or  density  functional  methods  to  the  MBS?  systems  and  using 
these  results  to  input  the  correct  level  and  type  of  Cl  required  for  the  system. 


Conclusions 

In  summary,  the  work  proposed  was  carried  out  with  a  large  degree  of 
success.  Excellent  new  parameters  have  been  developed  for  Si/  S,  and 

P  While  the  research  leading  to  SAMI  parameters  on  iron  are  not  comp  e  , 
L  is  onira  maS  of  tim/at  this  pofnt.  Substantial  obstacfos  have  b«n 
overcome  and  a  clear  route  to  the  objective  has  beeri  established.  The 
theoretical  framework  is  in  place  and  all  that  is  required  is  completion  of  the 

parameterization  process.  The  rate  at  which  such  but  ton 

Lried  out  was  perhaps  over-estimated  in  our  initial  but  iron 

should  be  finished  very  soon.  The  experience  we  have  gauged 
to  parameterize  SAMI  for  other  transition  metals  very  quickly..  The  primary 
quesrion  of  this  research  as  to  whether  the  Dewar-style  semiernpirical 
methodology  could  be  extended  to  transition  metals  has  been  answered  m  the 
very  definite  affirmative. 


Appendix  1:  SAMI  Results  for  Silicon 


Error 


Property 

#  Soecies 

SAMI 

AMI 

PM3 

AHf(kcal) 

62 

6.09 

7.33 

6.56 

IlI(D) 

20 

0.36 

0.41 

0.75 

IP<eV) 

43 

0.54 

0.58 

0.77 

Appendix  2:  SAMI  Results  for  Phosphorous 


Error 


Property 

#  Soecies 

SAMI 

AMI 

PM3 

AHf(kcal) 

60 

9.97 

16.63 

15.82 

11(D) 

17 

0.66 

0.99 

1.02 

IP  (eV) 

38 

0.69 

0.84 

0.93 

Appendix  3:  SAMI  Results  for  Sulfur 


Error 


Property 

#  Soecies 

SAMI 

AMI 

PM3 

AHf(kcal) 

76 

6.45 

7.09 

9.33 

\l(D) 

45 

0.40 

0.49 

0.59 

IP  (eV) 

67 

0.39 

0.33 

0.39 

Appendix  4:  Preliminary  SAMI  Results  for  Iron 

. ErrflrJm . •w . : 

■  m. . : - - Dipole  Moment 

-5.52  ^ 

-5.7 


FeH 

-3.11 

FeBra 

43.4 

FeBrs 

63.0 

FezCU 

10.4 

Fe2l4 

16.6 

FeCh 

3.2 

FeCb 

-16.2 

Fe(CO)4 

5.2 

Fe(CO)5 

-4.9 

Ferrocene 

-50.1 

FeCH3 

9.4 

FeOH 

-12.3 

Fe(OH)2 

34.6 

Average  (MU) 

30.5 

#  Molecules 

26 

-6.1 

-2.7 

-0.2 

2.2 

-6.3 


1.06  4.89 

15  2 


1.  kcal/mol.  2.  eV.  3.  Debyes. 


